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Abstract 

In order to generate the baryon asymmetry of the Universe sufficiently strong 
CP violation is needed. It was therefore proposed that at finite temperature 
there might be spontaneous (transitional) CP violation within the bubble walls 
at the electroweak phase transition in supersymmetric models. We investigate 
this question in the MSSM. 



1 Introduction 

For producing the baryon asymmetry of the Universe, we need extensions to the 
Standard Model. One of the Zakharov conditions requires nonequilibrium. In 
the MSSM this can befulfilled by a strong enough first order phase transition 
with a light scalar topU. Also a small bubble wall velocity-^ seems to support 
baryogenesis. But Zakharov's conditions also require CP violation and in the 
MSSM there are several mechanisms known to generate it. Explicit CP violating 
operators might conflict to experimental EDM boundJI. It were interesting to 
have a mechanism generating enough CP violation for baryogenesis without any 
conflict to experiments. While spontaneous CP violation is excluded at T = 
for the experimentally allowed parameter valued, there is a suggestion that it 
might be more easily realized at finite temperatureSH. 

Previously, the moduli of the two Higgs doublets around the phase boundary 
have been determined from the 2-loop effective potential!. The CP violating 
phase between thetwo Higgs doublets has been addressed perturbativelj&£9 and 
nonperturbativeljcJ. 

We present the first complete solution of the equations of motion for the 
phase between the two Higgs doublets within the MSSM, utilizing a perturbative 
effective potential, but without restricting it to the effective quartic couplings. 
Our conclusion^ differ from those obtained earlier on. 



2 Searching for CP violating phases 

We parameterize the two Higgs doublets of the MSSM as 

In addition, because of gauge invariance, the effective Higgs potential depends on 
the phases only via 9 = 9\ + 62, and we have an additional constraint h\d^6\ = 
h\d^82- We can then concentrate on 6. Assuming tree-level kinetic terms and 
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moving to a frame where the bubble wall is static and planar, the action to be 
minimized is 



r rl 1 1 h 2 h 2 i 

S oc J dz[-{d z h x f + -{d z h 2 f + --^-^{d z ef + F T (/»i, fca,0)] , (2) 

where Vr(hi,h2,9) is the finite temperature effective potential for hi,h2,8. In 
general, we are solving the equations of motion for hi,h2,9 following from this 
action. In the numerical solution we use the method outlined ini which deals 
with the minimization of a functional of the squared equations of motion. 

At the first stage, we consider the case with no explicit CP phases, and ask 
whether a particular solution without CP violation (9 = 0, it), is in fact a local 
minimum of the action or not. Clearly, it is not if 



2 / j i \ 1 d 2 V T (h 1 ,h 2 ,( 
m 3 {hi, h 2 ) - 



h,h 2 \ d9 2 



< 0, (3) 

8=0 



where we have divided by |/ii/i2|, assuming that this is non-zero. Eq. (|3|) is to 
be evaluated along the path found by solving the equations of motion for hi,h 2 - 
We have chosen the convention that h\ can have either sign, allowing us to 
consider only 8 = 0. For the case of the most general quartic two Higgs doublet 
potential, Eq. (||) agrees with the constraint on which most of the investigations 
of spontaneous CP violation are based. However, Eq. (Q) is true more generally, 
independent of the form of the potential Vr(hi, h,2,0). 
The tree-level potential of the theory is 

^tree = \m\h\ + -m\h\ + m^h^ cos9 + ^{g 2 + g' 2 ){h\ - h 2 2 ) 2 , (4) 

where g,g' are the SU(2) and U(l) gauge couplings, and at tree-level 

m i2 = -rm^ sin 2/3. (5) 

It follows that 771-3 ( /i i, /12) = (1/2)771-^ sin2/3 > 0, so that the minimum of the 
potential in the 9 direction is at 9 = 0. Thus, in order to get spontaneous CP 
violation one needs radiative corrections which can overcome the tree-level term. 

Older consideration]!!'! are based on the approximation to the effective po- 
tential where only the quadratic and quartic operators are considered. At finite 
temperatures around the electroweak phase transition, important contributions 
come from infrared sensitive non-analytic contributions which are not of this 
form, and can affect spontaneous CP violatior£J't3. Thus, it is important to 
solve the equations of motion more generally for the full effective potential. 

Here we consider the full finite temperature 1-loop effective potential of the 
MSSM. It is known that 2-loop corrections are very important in the MSSlvfl 
allowing for larger values of /ii,/i2 in the broken phase. Nevertheless, for the 
present problem we find that even 1-loop effects are in most cases very small, so 
we do not expect qualitative changes from the 2-loop effects. 



3 A scan for spontaneous transitional CP violation. 

The tree- level part of the effective potential Vr(hi, h 2 , 9) is in Eq. (Q). In the 
resummed 1-loop contribution to Vr(h\, h 2 , 9), we include gauge bosons, stops, 
charginos and neutralinos. This introduces dependences on the trilinear squark 
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mixing parameters At and n as well as on the squark mass parameters tjiq, mfj, 
and the U(l), SU(2) gaugino parameters Mi and M 2 . 

We now wish to see whether the constraint in Eq. (||) can be satisfied at the 
bubble wall between the symmetric and broken phases. To do so, we have to 
search for each parameter set for the critical temperature T c , solve the equations 
of motion for (hi, hi) between the minima, and evaluate m^(hi,hi) along this 
path. Since this is quite time-consuming, we proceed in two steps. 

1. At the first stage, we do not solve for hi,h2,T c , but rather take them as 
free parameters in the ranges hi/T = -2. .2 and h 2 /T = 0..2, T = 80. ..120 GeV. 
The zero temperature parameters are varied in the wide ranges 

tan/3 = 2. ..20, niA = 0...400 GeV, 

mu = -50.. .800 GeV, m Q = 50.. .800 GeV, (6) 
H,A t ,M 1 ,M 2 = -800.. .800 GeV. 

Here a negative mjj means in fact a negative right-handed atop mass parameter, 
— \mfj\. We have also studied separately the (dangerous^) region where the 
transition is very stronJ3'0, corresponding to mu ~ —70... — 50 GeV. 

Note that since we do not solve for the equations of motion at this stage 
but allow for h\ = ±\h\\, we have to divide in Eq. (||) by h\h 2 instead of 
\h\h 2 \: this leads in general to positive values due to the tree- level form of the 
potential, Eqs. (|]),(|]). A signal of a potentially promising region is then a small 
absolute value of the result, since this means that we are close to a point where 
dgVr(hi,h 2 ,9) crosses zero. 

2. At the second stage, we study the most favourable parameter region 
thus found in more detail. First of all, we search for the critical temperature. 
Then, we solve the equations of motion for (h\,h 2 ). By comparing with the 
exact numerical solution in several cases, we find that a sufficient accuracy can 
be obtained in practice by searching for the "ridge" as an approximation to 
the wall profile. It is determined as the line of maxima of the potential in the 
direction perpendicular to the straight line between the minima. Finally, we 
look for the minimum of Vr(hi, h 2 ,9) at fixed (h\, h 2 ): this is a fast and reliable 
approximation for the full solution in the case that 9 is small (i.e., just starts to 
deviate from zero), and corresponds to Eq. @. 

For the first stage, we perform a Monte Carlo scan with about 2 • 10 9 config- 
urations. Small values of m§(/ii, h 2 ) are scarce, and even then do not necessarily 
correspond to the desired phenomenon of spontaneous CP violation: they could 
also be points far from the actual wall. This can be clarified at stage 2. 

The parameter region found depends most strongly on rriA, tan/?, with a 
preference on small values of txia and large of tan/3, such that m\ 2 in Eq. $t|) 
is small. (This is in contrast to the requirements of a strong phase transition^) . 
There is also a relatively strong dependence on At and \i: the region favoured is 
shown in Fig. |]. The dependences on the other parameters are less significant; for 

re preferred. The region found is in rough agreement 



mu and mg small w 
with those found hit 

At the second stage, we make further restrictions. For instance, we exclude 
the cases leading to non-physical negative mass parameters. We also exclude 
cases leading to T = spontaneous CP violation in the broken phase: this 
phenomenon requires very small valued of m&- We also discard phase transitions 
which are exceedingly weak, v/T <C 0.1. 
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Figure 1 : The average value of m\ versus fj, and A t . We observe that small values of are not typical 
in any part of the plane but are on the average more likely for small fj,, Af, and that the distribution 
is wider (and thus more favourable) for like signs of fi,A t , as shown by the noisy contours obtained 
with a finite amount of statistics. 

ll the special point m% « was considered. Since il the thermal mass 
corrections were neglected this corresponds in the physical MSSM to a case 
where mfj + #T 2 ~ 0. Expanding the 1-loop cubic term from the stops to a 
finite order in v% jvi , it was suggested that transitional spontaneous CP violation 
can take place. This region isouite dangerous due to the vicinity of a charge 
and colour breaking minimunc3. Without expanding the 1-loop contribution in 
v\/v2, we cannot reproduce the behaviour proposed there. In any case, even 
before taking into account the experimental lower limits on the Higgs masses, 
we cannot find any promising case in the sample of ~ 2 x 10 6 configurations of 

StaE ; e 2 condudfl that after t aain g ftfto account the i„fta red sensitive effects 
inherent in the 1-loop effective potential, coming from a light stop and gauge 
bosons, and solving for the wall profile from the equations of motion, sponta- 
neous CP violation does not take place in the physical MSSM bubble wall. 
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